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Pulse trapping inside one-dimensional photonic crystal with relaxing cubic 

nonlinearity 

Denis V. Novitsk^ 
B.I. Stepanov Institute of Physics, National Academy of Sciences of Belarus, 
Nezavisimosti Avenue 68, 220072 Minsk, Belarus. 
(Dated: May 12, 2010) 

We theoretically study the effect of pulse trapping inside one-dimensional photonic crystal with 
relaxing cubic nonlinearity. We analyze dependence of light localization on pulse intensity and 
explain its physical mechanism as connected with formation of dynamical nonlinear cavity inside 
the structure. We search for the range of optimal values of parameters (relaxation time and pulse 
duration) and show that pulse trapping can be observed only for positive nonlinearity coefficients. 
We suppose that this effect can be useful for realization of optical memory and limiting. 
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I. INTRODUCTION 

Photonic-band-gap structures are actively studied as 
promising elements for different devices of nonlinear and 
quantum optics [l|-|3| . Study of nonlinear photonic crys- 
tals is connected with the possibility of dynamical adjust- 
ment of the parameters of the system. Photonic crystals 
(including two- and three-dimensional) allow to control 
dispersion and diffraction properties of light [J, [a], obtain 
pulse reshaping 6], pulsed high-harmonic generation [7J, 
etc. The processes of light self-action in such systems re- 
sult in localization effects such as gap solitons formation, 
discrete mode existence or light energy localization (see, 
for example, |a-lllj). 

The simplest example of one-dimensional photonic 
crystal can be represented as a periodic set of alternating 
dielectric layers with large depth of refractive index mod- 
ulation. One of the most prominent effects of nonlinear 
optics of such structures is strong pulse compression in 
photonic crystals with non-resonant cubic nonlinearity 
[13 . Il3j . The effect of light pulse localization, or trap- 
ping, in one-dimensional photonic crystal with a defect 
was studied in Ref. |14| and, in more general form, in 
Refs. [H[l3. 

In this paper we study pulse propagation in photonic 
crystal with relaxing cubic nonlinearity. It is obvious that 
if we reduce incident pulse duration, the inertial proper- 
ties of medium nonlinearity should be taken into account. 
Indeed, the lowest values of relaxation time connected 
with the electronic Kerr mechanism (a few femtoseconds 
17]) appear to be comparable with pulse durations ob- 
tained at the modern set-ups. As it was shown in Ref. 
[l8| , relaxation of nonlinearity results in vanishing of the 
effect of femtosecond pulse compression which can be ob- 
tained in relaxation-free case. 

In our research we use numerical simulations of the 
Maxwell wave equation taking into account the process 



of nonlinearity relaxation. The method used allows to 
obtain numerical solutions of the problem without any 
assumptions about medium parameters modulation or 
rate of variation of field envelope. We show that, in a 
certain region of pulse amplitudes, relaxation times and 
pulse durations, the pulse can be trapped inside a nonlin- 
ear cavity dynamically formed by the light. Appearance 
of the cavity is connected with local change of reflective 
properties of the photonic structure. In other words, dy- 
namical shift of band spectrum occurs. Trapping pulse 
inside the photonic crystal is a prospective effect for such 
possible applications as optical limiters (which transmit 
light only with proper intensity), optical buffers or mem- 
ory (which allows to store light for time, large in com- 
parison with characteristic transmission time). Large ad- 
vantage of the scheme considered is absence of necessity 
to introduce nonuniformity or any imperfections in the 
structure of nonlinear photonic crystal. 

The article is divided into several sections. In Section 
[IT] the main equations are given and the approach for nu- 
merical solving the wave equation is considered. Section 
IIIII is devoted to some phenomenological aspects of pulse 
trapping effect. In Section IIVI the physical mechanism 
of pulse trapping in photonic crystal with relaxing cu- 
bic nonlinearity is discussed. Finally, in Section [V] we 
consider some conditions for pulse trapping observation 
connected with the proper choice of pulse duration and 
relaxation time. 



II. MAIN EQUATIONS AND NUMERICAL 
METHOD 



In this paper we consider ultrashort pulse interaction 
with one-dimensional photonic crystal made of substance 
with relaxing cubic (Kerr) nonlinearity. Light propaga- 
tion along z-axis is governed by the Maxwell wave equa- 
tion 
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where E is electric field strength, n is medium refractive 
index that depends on light intensity I = \E\ 2 as 



n = n Q (z) + 5n(I, t). 



(2) 



Here tiq(z) is a linear part of refractive index. Time de- 
pendence of nonlinear term Sn is responsible for relax- 
ation process and is described by the first-order differ- 
ential equation due to the Debye model of nonlinearity 
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where n 2 is Kerr nonlinear coefficient, t n i is the charac- 
teristic relaxation time. We consider fast relaxing me- 
dia (electronic Kerr mechanism) with relaxation times as 
small as few femtoseconds. Representing field strength 
as E = A(t, z) exp [i{ujt — kz)], where w is a carrier fre- 
quency, k — lj/c is the wavenumber, and introducing 
new, dimensionless arguments t = tut and ^ — kz, 
we come to the wave equation for the pulse amplitude 
A(t,z), 
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+ (n 2 -1)A = 0. (4) 



Usually, second-order derivatives are neglected at this 
point resulting in slowly varying envelope approximation. 
However, in the case of abrupt changes of refractive index 
(photonic crystal) it may fail, so that we should solve the 
full Eq. (01 . In addition, this equation allows to describe 
the behavior of electric field in the structure without divi- 
sion into forward and backward waves. Equation Q and 
the computational scheme considered below are similar 
to those of Refs. [20, [21( where they were implemented 
to consider light propagation in a dense resonant medium 
and a photonic crystal containing it. 

Equation (Q| can be solved numerically by using the 
finite-difference time-domain (FDTD) approach. The 
computational scheme is based on calculation of ampli- 
tude value at every mesh point (ZAr, j'A£) as 

4 +1 - [-a x A 1 - 1 + b x A l j+x + b 2 A\_ l + fA]]/a 2 . (5) 

Here the auxiliary values are 

oi = (n l f l ) 2 (l - jAr), a 2 = (n l + l f{l + iAr), 
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/ = 2(^.) 2 -2^J +Ar 2 [(4) 2 -1]. 

The values of refraction index at the mesh points can 
be obtained in terms of finite-difference representation of 
Eqs. ©-©, 



;/ , ' = noj + 5n l+1 , 
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FIG. 1: Scheme of a photonic crystal considered. Parameters 
of it: refractive indices n a = 2, rib = 1.5; thicknesses a = 0.4, 
b = 0.24 /im; number of layers N — 200. 
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FIG. 2: (Color online) Interaction of ultrashort pulse with 
the photonic crystal with (a) relaxation-free (t n l = 0) and 
(b) relaxing (t n i = 6 fs) cubic nonlinearity. Pulse peak in- 
tensity I m = -^o, where 7o is such that 71.2/0 = 0.005. Other 
parameters are the same as in the caption of Fig [TJ and in the 
text of the paper. The inset demonstrates the spectra of the 
photonic crystal (solid line) and of the pulse (dotted line). 



where r„; = u)t n i , hqj and n 2 j are the values of the back- 
ground refraction index and nonlinearity coefficient at 

To correctly set the boundary conditions we use the 
total field / scattered field (TF/SF) method and the per- 
fectly matched layer (PML) method which allows to ap- 
ply the so-called absorbing boundary conditions at the 
edges of the calculation region [lj| . 



III. PULSE TRAPPING EFFECT 

Let us consider propagation of ultrashort (femtosec- 
ond) light pulse in one-dimensional photonic crystal 
shown in Fig. [T] Spatial periodic modulation of the 
background refractive index no (z) defines the structure 
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FIG. 3: (Color online) (a) Transmitted and (b) reflected in- 
tensity after interaction with nonlinear photonic crystal at 
different incident pulse peak intensities (A m — 1, 3, 7 in units 
of Ao). Relaxation time of nonlinearity t n i = 6 fs. Other 
parameters are the same as in the caption of Fig [2] 
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FIG. 5: Distribution of light intensity inside the photonic 
crystal at different time points. Pulse peak amplitude (a) 
A m = 3Ao, (b) Am = 7Aq. Relaxation time of nonlinearity 
t n i — 6 fs. 
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FIG. 4: (Color online) Dependence of transmitted, reflected 
and overall output energy (as fraction of input energy) on 
incident pulse peak amplitude A m . Relaxation time of non- 
linearity t n i — 6 fs. Energy was integrated over the time 
interval of 200£„. 



of it, so that it can be treated as a set of alternate lay- 
ers. The parameters used in calculations are as follows: 
refractive indices of the layers n a = 2, nj, = 1.5; their 
thicknesses a = 0.4, b = 0.24 /im; number of layers 
TV = 200. Nonlinear coefficient of the material is defined 
as 712-Z0 — 0.005, i.e. the pulse amplitude is normalized 
by the value Aq = \/7o- This value of cubic coefficient 
provides refraction index change of thousandth and hun- 
dredth of unity and was used in calculations of Ref . [18| . 
Note that nonlinear coefficient is positive (focusing non- 
linearity), so that refractive index increases with the in- 
tensity. The incident pulse is assumed to have Gauss 
envelope A — A m exp {—t 2 /2ti). Here t p is a pulse dura- 



tion which further takes on the value of 30 fs, while the 
carrier frequency lies on the wavelength A = 1.064 fim. 

Figure [2] shows the results of calculations of pulse inter- 
action with nonlinear photonic structure with and with- 
out relaxation. One can easily see that the effect of pulse 
compression (obtained at t n i = 0) is completely absent 
when t n i = 6 fs. This result is in strict accordance with 
the conclusions of Ref. [18|] . 

Vanishing of compression effect in Fig. [2] was ob- 
tained for pulse peak intensity I m = \A m \ 2 = I , where 
Iq = \Aq\ 2 is the value of intensity corresponding to 
ri2lo = 0.005. If we take greater A m , the transmitted 
pulse continues to decrease. As it is seen in Fig. |3ja), 
only small part of incident light can pass through the non- 
linear photonic crystal at A m = ZAq. This situation is 
observed at A m = 7A as well, though the reflected pulse 
gets larger in this last case [Fig. [3jb)]. When we inte- 
grate intensity of reflected and transmitted light over a 
certain time interval (200i p in our calculations, i.e. large 
enough in comparison with time required for pulse to pass 
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FIG. 6: Distribution of light intensity inside the photonic 
crystal at different time points. Pulse peak amplitude A m — 
5Aq. Relaxation time of nonlinearity t n l = 0. 



through the photonic crystal which is about 30£ p , see Fig. 
[5]), we obtain the characteristic energy curves (Fig. [3} 
that describe the change in behavior of pulses as their 
peak intensity is increasing. It is seen that, when A m 
gets larger than 2A$ : reflected, transmitted and overall 
(summarized) energies demonstrate dramatic decrease. 
After reaching the minimum (output energy is about 20 
percent of the input one), the curves begin to rise slowly. 
This slow growing of transmitted energy continues, while, 
for reflected and total ones, rise becomes more steep (es- 
pecially, when A m > 5Aq) and, finally, the plateau is 
observed for A m > 7Aq. Thus, the input energy almost 
entirely transforms to reflected light at high intensities of 
incident pulse. 

It is obvious that the energy of pulses with amplitudes 
between A m « 2Ao and A m ss 6Aq is confined inside 
the photonic crystal. Figure [5][a) demonstrates light in- 
tensity distributions along the length of the structure at 
different instants of time for the pulse peak amplitude 
A m — 3Aq. At first (t — 50£ p ), the most part of pulse 
energy is localized in a narrow region of the nonlinear 
photonic crystal, near the position L = 50 ^m (the total 
length of the structure is about 130 /im). As time goes 
by, energy tends to redistribute more uniformly. Never- 
theless, there is still pronounced maximum of intensity 
distribution, moreover it is shifted towards larger posi- 
tions, namely 60 — 70 /im at t > 700t p . At large time 
points this distribution stays almost invariant, or station- 
ary. Its maximum only slightly decreases, which seems to 
be connected rather with further redistribution than with 
output radiation. Anyway, even at t = 4000£ p approxi- 
mately 80 percent of pulse energy is still confined inside 
the photonic structure, just as at t = 200t p (see Fig. H|). 
This time is more than 100 times higher than the interval 
needed for pulse to pass through the system. Recalling 
that t p — 30/s, this delay time in absolute units is greater 
than 100 ps. Therefore we can say about pulse trapping 



in this case. Only for t > 4000£ p the system starts to 
slowly emit light so that the sharp distribution shown in 
Fig. [SJa) becomes violated. More detailed calculations 
show that even at t = 10000i p about 50 percent of the 
initial energy is still inside the photonic crystal, though 
it is distributed much more uniformly. 

Now let us consider the pulse with A m = 7A$ . The cor- 
responding spatial distributions are shown in Fig. [5jb). 
It turned out that in this case the pulse is localized near 
the very beginning of the structure. In this position it 
rapidly loses energy which is mainly radiated through 
the front (input) end of the system. So this radiation 
gives significant contribution to reflection [see Fig. [3jb)]. 
After light intensity becomes lower than a certain thresh- 
old, the pulse starts moving and widening. It moves quite 
slowly, so that some part (about 20 percent) of the input 
energy is confined inside the photonic crystal for a long 
time, but the peak intensity of the distribution is very 
low if we compare with the case of A m = 3Aq. 

It is interesting to compare pulse behavior consid- 
ered with pulse propagation in photonic crystal with 
relaxation-free nonlinearity, i.e. at t n i = 0. Intensity 
distributions inside the system for this case are demon- 
strated in Fig. |6l It is seen that there is no any long-term 
energy localization in non-relaxing photonic structure. 
Light exhibits only chaotic "wander" inside it and simul- 
taneous attenuation due to emitting through input and 
output ends. Finally, almost all energy of the pulse is 
already radiated by the instant of time t = 700i p . There- 
fore we can say that relaxation of nonlinearity is a nec- 
essary condition to obtain the effect of pulse trapping 
inside a photonic crystal. 



IV. PHYSICAL MECHANISM OF PULSE 
TRAPPING 

What is the physical reason, or mechanism, of this phe- 
nomenon? As pulse propagates inside the photonic crys- 
tal possessing cubic nonlinearity, refractive index of the 
structure changes dynamically according to light inten- 
sity. If nonlinearity is relaxation- free, these changes are 
instantaneous and depend entirely on field distribution 
at current instant of time. In the case of relaxing non- 
linearity, nonlinear variation of refractive index (Sn in 
our notation) can form a certain stable structure due to 
retardation in its change. Appearance of this nonlinear 
dynamical " cavity" results in pulse trapping: light tends 
to change the distribution of Sn and leave the cavity, but 
inertia of nonlinearity stabilizes it so that intensity is 
transformed to provide steady spatial distribution 8n(z). 
The example of this is shown in Fig. EJb). One can 
see that maximal values of nonlinear variation of refrac- 
tive index Sn (so-called primary maxima) are achieved in 
the layers with low linear refractive index uq = 1.5 [see 
Fig. EJa)]. This is in accordance with the effect of light 
concentration in low refractive index regions of photonic 
crystal for radiation tuned to high-frequency side of re- 
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FIG. 7: Spatial variation of (a) linear part of refractive index 
rio(z), (b) nonlinear contribution Sn. Pulse peak amplitude 
A m — 3Ao, time point t = 50£ p . Relaxation time of nonlin- 
earity t n i — 6 fs. 
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FIG. 9: Spatial variation of (a) intensity, (b) nonlinear vari- 
ation 5n in the case of uniform cubic medium with linear 
part of refractive index no = 1.8125. Pulse peak amplitude 
A m — 20^4.0, time point t — 1000£ p . Relaxation time of non- 
linearity tnl = 6 fs. 
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FIG. 8: (Color online) Transmitted and reflected pulses after 
interaction with a uniform cubic medium. Pulse peak ampli- 
tude (a) A m = 3^o, (b) Am — 20 Aq. Linear part of refractive 
index no = 1.8125. Relaxation time of nonlinearity t n l = 6 
fs. 



flection spectrum |2j as it is in the case considered (see 
the inset in Fig. [2]). On the other hand, high refractive 
index layers of the structure contain minor (secondary) 
maxima of Sn. 

Obviously, low-intensity initial pulses are insufficient to 
produce high enough nonlinear contribution Sn, so that 
pulse is mainly transmitted and reflected during short 
time after pulse incidence. For higher intensities [e.g. 
A m = 3Aq as in Fig. Ufa)] pulse forms nonlinear cav- 
ity inside the photonic crystal which localizes the most 
part of pulse energy. Though this cavity slowly tends to 
uniformity of Sn and loses its energy, it allows to trap 
pulse light for relatively large time interval. This behav- 



ior is obtained in wide range of pulse amplitudes, approx- 
imately from 2.5Aq to 5Aq. As pulse intensity increases, 
cavity formation position moves towards the front end 
of the photonic structure. Finally, for A m « 7Aq the 
nonlinear cavity appears so close to the entrance of the 
system [Fig. Hfb)] that it rapidly emits almost all its 
energy in the form of reflected light. That is why only 
intermediate intensities of pulse (not too low and not too 
high) are suitable to obtain the effect of trapping pulse 
inside the structure considered. 

Finally, we should explain the role of photonic crystal 
in this process. Is it necessary to use photonic-band-gap 
structure or, maybe, the effect of pulse trapping can be 
observed in a uniform cubic medium? To clarify this 
question we performed calculations of pulse interaction 
with the medium possessing relaxing nonlinearity and 
mean refractive index rio(z) = (an a + brib)/(a + b) = 
1.8125. One can see [Fig. [Ha)] that in this case there is 
no any light localization for pulse amplitude A m = 3Ao. 
Only much more intensive pulses start to lose consid- 
erable part of energy inside the medium. For example, 
about 40 percent of pulse energy is found to be trapped 
for A m = 20 A [Fig. [Sfb)]. This part of energy stays in- 
side the uniform medium for a long time and is connected 
with corresponding nonlinear variation of refractive index 
(see Fig. [9] for the instant of time t = 1000t p ). However, 
the distribution of light intensity in the uniform medium 
seems to be stochastic and does not resemble pulse en- 
velope as in the case of nonlinear photonic crystal [com- 
pare the distributions in Figs. Ufa) and[5ja)]. Therefore 
we cannot call light localization in the uniform medium 
with relaxing nonlinearity by the pulse trapping in the 
full sense of this term. Moreover, in photonic crystals 
we need pulse intensities which are less by an order of 
magnitude than in the case of the uniform medium. 

Perhaps, the nonlinear cavity formation is connected 
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FIG. 10: (a) Distribution of nonlinear refractive index vari- 
ation 8n corresponding to intensity distribution of Fig. [5ja) 
at t = 1000tp. (b) Difference between reflection coefficients 
of the partial structures including the layers from the input 
to a certain final position in the cases of Sn given by picture 
(a) and Sn — 0. Calculations were carried out for the central 
wavelength A = 1.064 /itn. 



with local change of reflective properties of photonic crys- 
tal. This results in dynamical shift of the band spec- 
trum of the structure. To examine this situation we 
consider the refractive properties (for the central wave- 
length A = 1.064 /im) of the structure with refractive 
index variations shown in Fig. [Tola) [it corresponds to 
intensity distribution of Fig. EJa) at t = 1000i p ]. We 
calculate reflectivity of the partial structures, which in- 
clude the layers from the input to a certain final position 
(inside the whole photonic crystal) , and compare it with 
the case of linear structure. Difference between reflectiv- 
ities in this two cases as a function of final position is 
demonstrated in Fig. [TUJb). It is seen that the struc- 
ture with modified refractive index modulation provides 
locally large reflectivity deviations from the linear case. 
This deviations mainly appear at large positions inside 
the crystal where 5n is high enough, so that they can 
prevent light propagation in forward direction. Similar 
situation is observed if we consider backward propaga- 
tion. Hence, it turns out that light is trapped in the 
central region of the structure. Finally, it is worth not- 
ing that this picture of refractive index variations and 
reflectivity deviations is permanently changing, though 
it is stabilized by the relaxing properties of nonlinearity. 



V. ON OPTIMAL CONDITIONS OF PULSE 
TRAPPING 

In previous sections we considered the effect of pulse 
trapping in photonic crystal with relaxing cubic nonlin- 
earity only for a single set of time parameters: t n i = 6 
fs and t p = 30 fs. Fig. ITTT a) shows the dependence of 
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FIG. 11: (Color online) Dependence of transmitted, reflected 
and overall output energy (as fraction of input energy) on (a) 
relaxation time at t p = 30 fs, (b) pulse duration at t n i = 10 
fs. Incident pulse peak amplitude A,„ — 4Ao- Energy was 
integrated over the time interval of 200£ p o , where t p o — 30 fs. 
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FIG. 12: Dependence of overall output energy on the ratio 
tnl/tp plotted according to the data of Fig. QT] 



output energies (reflected, transmitted and total) on re- 
laxation time at fixed value of pulse duration t p — 30 fs. 
The behavior of these dependencies is similar to that of 
the curves in Fig. |4j abrupt drop in the range of small 
t n i (less than 1 fs), smooth increasing and, finally, steep 
rise of the curves for reflected and total energies. The 
full range of relaxation times where the pulse trapping 
can be observed is rather wide: from a fraction of fem- 
tosecond (relaxation is so fast that it does not influence 
the pulse) to about 150 fs that is much greater than t p 
(medium reacts so slowly that the nonlinear cavity forms 
near the very entrance of the system) . The optimal value 
is tni ~ 10 fs. 

If we fix relaxation time t n i — 10 fs and vary pulse du- 



1.0 



0.8 



03 r, a 

c 0.6 



a> 0.4 

N 

"re 

E 

O 0.2 



0.0 



incident 



reflected 



transmitted 



— r* 
25 50 

Timef/f 

p 



75 



much more wide range than pulse duration. This also im- 
plies that we can use media with relatively slow relaxing 
nonlinearities to obtain the effect of pulse trapping. 

Another question is connected with the role of sign 
of nonlinearity. So far we considered only positive non- 
linearity coefficient such that 712/0 = 0.005. If we take 
772/0 = —0.005 (defocusing nonlinearity), there is no any 
symptoms of pulse trapping as one can see in Fig. 1 131 for 
the amplitude A m — 3Aq. Though this problem should 
be studied in detail, the preliminary conclusion is that 
trapping can be observed only for 712 > 0, at least for 
comparatively low intensities. 



VI. CONCLUSION 



FIG. 13: (Color online) Transmitted and reflected pulses after 
interaction with a photonic crystal with negative nonlinearity 
coefficient 712/0 = —0.005. Pulse peak amplitude A m = "AAq. 
Relaxation time of nonlinearity t„i = 6 fs. Other parameters 
are the same as in the caption of Fig fT] 



ration, the behavior of output energies is quite different 
[Fig. ITTTb)]. As in the previous case, it demonstrates 
abrupt decrease for small times, namely, pulse durations 
t p ~ 10 fs which correspond to only a few optical cy- 
cles. Such very short pulse is not able to create a stable 
nonlinear cavity. On the other hand, if t p is increasing, 
there is only stepless growth of output energies, reflected 
and transmitted ones being approximately of the same 
magnitude. The reason is the same: nonlinear cavity is 
not formed as long as a whole pulse cannot be placed 
inside the structure. The optimal value of pulse duration 
is about 20 fs. 

Since the limits of pulse trapping region are due to dif- 
ferent reasons in the cases of fixed t p and t n i , the width 
of this region will be different, too. In Fig. Q2]we plotted 
the curves of Fig. [TT] versus the ratio t n i/t p . This fig- 
ure shows explicitly that relaxation time can be varied in 



To summarize, we have analyzed the possibility of 
trapping pulse in photonic crystal with relaxing cubic 
nonlinearity. By using numerical simulations, we showed 
that this process is due to the balance between light 
spreading and inertia of nonlinearity which results in 
steady nonlinear cavity formation and pulse trapping 
within it. Photonic crystal is a necessary element for this 
cavity to appear, due to the processes of dynamical lo- 
cal change of reflection and transmission of the nonlinear 
structure, and, in addition, leads to decreasing of pulse 
intensity required to observe trapping. We discussed the 
reasons for pulse trapping disappearance at high and low 
values of both pulse duration and relaxation time result- 
ing in existence of the range of optimal magnitudes of 
these parameters. 
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